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Abstract. We report proper motion measurements for 25 
very-low mass (VLM) star and brown dwarf (BD) candi- 
dates of the Pleiades cluster previously identified by Bou- 
vier et al. (1998). Proper motions are measured with an 
accuracy of 9 mas/yr, compared to an expected tangen- 
tial motion of about 50 mas/yr for Pleiades members. Of 
the 25 candidates, 15 have a membership probability of 
95% or more and 7 are rejected as being field dwarfs. The 
3 remaining candidates exhibit independent evidence for 
membership (lithium absorption or long-term proper mo- 
tion). From the firm identification of Pleiades VLM and 
BD members, the cluster's substellar mass function is re- 
vised to dN/dM (x M" 5 in the mass range from 0.04 to 

0. 3M Q . 

Key words: Stars: low-mass, brown dwarfs - Stars: lumi- 
nosity function, mass function - (Galaxy:) open clusters 
and associations: individual: Pleiades 

1. Introduction 

The determination of the stellar mass function is an im- 
portant challenge for various domains of astrophysics such 
as, e.g., the star formation process, the structure and evo- 
lution of the Galaxy, the dynamical evolution of stellar 
clusters and stellar systems, etc... This function describes 
the relative number of stars par unit mass and is usually 
approximated by power-law segments of the form dN/dM 
oc M _Q in various mass domains. While its shape is rel- 
atively well constrained from solar-type to massive stars 
{a ~ 2.3-2.7 for M>1 M , Salpeter 1955, Scalo 1998 and 
references therein), it is more uncertain for low-mass stars 
(a ~ 1.0-2.0 for 0.3<M<1 M , Kroupa 2000 and refer- 
ences therein) and, up to a couple of years ago, was still 
unexplored at very low masses in both the stellar (0.08<M 
<0.3 M ) and substellar (M<0.08 M ) regimes. The de- 
termination of the substellar mass function is indeed one 
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of the main motivations for the recent and exploding quest 
for free-floating brown dwarfs (see, e.g., Oppenheimer et 
al. 2000 for a review). 

An estimate of the mass function (MF) in the upper 
part of the substellar domain (~0. 04-0. 08 M ) now exist 
for the Pleiades cluster (Bouvier et al. 1998: a ~ 0.6±0.15, 
Martin et al. 1998, 2000: a ~ 0.5-1.0, Hambly et al. 1999: 
a ~ 0.7) and has recently been derived for field brown 
dwarfs in the solar neighbourhood (Reid et al. 1999: a ~ 
1.0-2.0 depending on the local stellar birthrate). Deter- 
minations of the substellar mass function have also been 
reported for several star forming regions (e.g. Comeron et 
al. 2000, Luhman & Rieke 1999, Hillenbrand & Carpenter 
2000) with, however, the additional difficulty of varying 
extinction within the molecular cloud which makes more 
problematic the derivation of an unbiased sample of very 
low-mass young objects that can be used to build the sub- 
stellar MF. 

Estimates of the substellar mass function of the 
Pleiades cluster rely on deep, wide-field photometric sur- 
veys that identify substellar candidates from their location 
in optical color-magnitude diagrams (CMD). One of the 
major concern regarding these samples of very low-mass 
(VLM) and brown dwarf (BD) candidates is the degree of 
contamination by foreground late-type field dwarfs, which 
may lie in the same region of the CMD as the Pleiades 
low-mass members. Bouvier et al. (1998) thus estimated 
at about 30% the level of contamination of their candi- 
date sample by field dwarfs. Some of the BD candidates 
have actually been confirmed through the "lithium test" 
(Rebolo et al. 1992) but, even with the largest telescopes, 
only the brightest substellar candidates are amenable to 
this test (Stauffer et al. 1998). Other Pleiades BD can- 
didates originally identified from their location in an (I, 
R-I) CMD have later been rejected on the basis of their 
discrepant near-IR colors (Martin et al. 2000). However, it 
is currently unclear whether all contaminating field stars 
can be unambiguously recognized in a near-IR CMD. 

A powerful way to recognize true Pleiades members 
among the photometric candidates is to measure their 
proper motion. Proper motion studies should allow us 
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to pick out members with a high degree of confidence 
because the cluster's peculiar motion is large compared 
to the non-member field stars (/i Q cos<5 — 19.15 ± 0.23 
mas/yr, (j, 5 = -45.72 ± 0.18, Robichon et al. 1999) and 
the intrinsic velocity dispersion of the cluster is small (~ 1 
mas/yr, Van Leeuwen 1980). Such a large motion can be 
measured from sets of sharp images separated by only a 
few years and as soon as a candidate is recognized as a 
Pleiades member on the basis of its kinematics, its status 
(VLM star or BD) directly follows from its photometric 
properties. 

We therefore obtained in September 1999 new images 
for 25 of the 26 VLM star and BD candidates of the 
Pleiades originally observed by Bouvier et al. (1998) in 
December 1996. In Section 2, we describe how proper mo- 
tion was derived for each candidate from the 2 sets of im- 
ages separated by nearly 3 years. In most cases, the results 
presented in Section 3 allow us to unambiguously identi- 
fied Pleiades members and they are compared with other 
diagnostics of Pleiades membership. In Section 4, we dis- 
cuss whether previous estimates of the Pleiades substellar 
function has to be revised in the light of these new re- 
sults and briefly discuss the kinematics of very low-mass 
members of the cluster with respect to their formation 
mechanism. 

2. Observations and astrometric reduction 

The first set of images, from which VLM and BD can- 
didates were identified, was obtained by Bouvier et al. 
(1998) in December 1996 with the CFHT 8K wide-field 
camera. Of the 26 identified candidates, 25 were reob- 
served in September 1999 with the CFHT 12K camera. 
Both instruments have the same pixel size (0.21") and 
the images were obtained under comparable seeing condi- 
tions (0.9"and 0.7", respectively). In 1999, the exposure 
time was adjusted so as to provide similar signal-to-noise 
ratio on the candidates as on the 1996 images. Astrometry 
was performed on I-band images from the 2 epochs. 

The principle of the astrometric procedure we used to 
measure the proper motion of the candidates is as fol- 
lows. We first measure the position (xg6,y96), (xg^ygg) 
of the candidate on the two sets of images using 
IRAF/CENTER. On the same images, we also measure 
the positions of typically 10 point-like objects, presum- 
ably background field stars, located within an angular ra- 
dius of about 3' from the candidate on the CCD. These 
stars are then used as relative astrometric references to 
compute the spatial transformation function that maps 
1999 coordinates to the 1996 reference frame using the 
IRAF/IMMATCH package. The rms uncertainty associ- 
ated to the spatial transformation, a t , is computed by the 
IRAF task as well as the residuals for each astrometric 
reference star. In case the residuals for a reference star is 
significantly larger than a t , the object is discarded and 
the transformation recomputed. 



This transformation is applied to (xgg,ygg) in order 
to project the 1999 coordinates of the candidate in the 
1996 reference frame: (x 99 ,y 99 ) — > (x 99 _g 6 , y 9 g^g 6 ). The 
displacement of the object in pixels between 1996 and 1999 
is simply given by: 

Ax = X 9 9^ 96 - X 96 = (At)/I x 
Ay = J/99^96 - 2/96 = (At)/i y 

where At is the time lag between the two epochs and 
the relative proper motions along the x and y 
axes of the CCD in the 1996 reference frame. The precision 
of the measurement is o\ x = <j\ y = 2<r^ y + a^, where a x . y 
is the rms error on the measurement of the position of the 
candidate on the CCD and a t the error associated with 
the spatial transformation from one epoch to the other. 
We evaluated a x . y ~ 0.07 pixels by measuring the posi- 
tions of the same objects on consecutive I-band images of 
the same field. Tests were run with the IRAF/GEOMAP 
package in order to minimize at- We found that a third 
degree polynomial transformation yields the best results 
with a t — 0.05 pixels. We verified a posteriori that the 
proper motion measured for each candidate was not undu- 
ely sensitive either to the degree of the polynomial trans- 
formation or to the number of astrometric reference stars 
used to compute the transformation. 

The displacement (/i x ,/i y ) in pixels/yr has still to be 
converted into a displacement on the sky (fj, a cosS, 
in mas/yr by calibrating the spatial scale and orienta- 
tion of the 1996 coordinates system. The calibration was 
done using PRIAM (Procedure de Reduction d'Images As- 
troMetriques), a software package developed by A.Fienga 
and J.Berthier at the IMCCE (Institut de Mecanique 
Celeste et de Calcul des Ephemerides). This software cal- 
culates the transformation which converts pixel positions 
(x, y) into J2000 celestial coordinates (a, 8) for a set of 
astrometric standards. For each candidate field, we mea- 
sured on the image the position of typically 10 stars be- 
longing to the USN02 catalogue whose astrometric ac- 
curacy is of order 500 mas. PRIAM then computes the 
coefficients A0,..,B2 of the transformation: 

a = A0 + A1 x + A2 y 
5 = B0 + Bl x + B2 y 
from which we obtain: 

\i a cos S = Al \i x + A2 fiy 

Us = Bl fi x + B2 fiy 

The rms error on these coefficients is ~ 0.1 mas/pixel, 
and can be neglected compared to a x . y and a t - Hence, the 
final uncertainty on the proper motion measurement is 

2 2 p%X . 2 2 \ 

oleosa ~ °w ~ TwJ y la x,v + a t ) 
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which typically amounts to 9 mas/yr rms for a pixel 
scale pix = 205 mas and At — 2.8 yr. 

3. Results 




fi a cos6 mas/yr 

Fig. 1. Proper motion vector point diagram for objects 
located in Pleiades fields. Filled dots show VLM and BD 
candidates and open circles the field stars used as astro- 
metric references. North is up and East is left. The cluster 
proper motion is indicated by a cross and two circles cen- 
tered on this value are drawn. The radius is equal to 2a 
and 3cr for the small and the big circle respectively where 
a corresponds to the 2D-gaussian dispersion of the candi- 
dates in this diagram (see text). 

The proper motions of the 25 Pleiades VLM star and 
BD candidates are listed in Table [l] and illustrated in a 
vector point diagram (VPD) in Figure |l|. The field stars 
used as astrometric references scatter around the origin of 
the diagram since they have been selected on the basis of a 
negligible proper motion. In contrast, most of the Pleiades 
candidates fall within 3u of the expected cluster's mean 
motion. Following the method outlined by Sanders (1971) 
we compute a membership probability for each candidate 
as follows: 

f c (Ha cos S, (is) 



ff(Ha COS S, (1$) + fc(fJ>a COS 5, fXs) 

where f c (p a cos 5, /is) and ff(/i a cos 5, fis) are the vector 
point distributions of the members and field stars respec- 
tively. We assumed f c is a bivariate normal function and 



we found 
f c (fi a cos6, us) = 0.04exp| - - ( 

'/i* + 40.1x2 



Ha cosS — 21. 5\ 2 
7J 



8.5 



01} 



by fitting a 2-D gaussian, plus an assumed uniform dis- 
tribution for field stars, on the vector-points distribution 
of the VLM and BD candidates. The field stars distribu- 
tion is then fj — 3.31CP 4 . We also tried to reproduce the 
distribution of both Pleiades candidates and astrometric 
references in the VPD by adding a 2-D gaussian centered 
on (0,0) to the fitting function. This leads to a dispersion 
of ~ 4.5 mas/yr for the astrometric references. Member- 
ship probabilities computed for the candidates by either 
method are very similar because the distributions of can- 
didates and astrometric references do not overlap in the 
VPD. 

The membership probability of the candidates is listed 
in Table 0. We also list in this table other diagnostics of 
Pleiades membership obtained by Stauffer et al. (1998) 
and Martin et al. (2000), namely: EW(Li), EW(H„), ra- 
dial velocity, and (I-K) index. We consider the 15 candi- 
dates for which we found p > 95% as Pleiades members. It 
is interesting to note that two of these confirmed members 
(CFHT-PL-7 a very low-mass Pleiades star, and CFHT- 
PL-25, the least massive brown dwarf) do not exhibit Ha 
in emission according to Martin et al. (2000). In addition 
to these 15 highly probable Pleiades members, CFHT-PL- 
12 has a membership probability of 88% and lithium has 
been detected in its spectrum (Stauffer et al. 1998) which 
makes it a bona fide Pleiades brown dwarf. We could not 
measure the proper motion of CFHT-PL-3 (= HHJ 22) 
because it is located both at the edge of a CCD and at 
the edge of the camera field of view which makes the com- 
putation of the spatial transformation between the two 
epochs unreliable. However, the long-term proper motion 
of HHJ 22 has been measured by Hambly et al. (1993) with 
an epoch difference of nearly 40 years and indicates highly 
probable Pleiades membership. We therefore consider this 
candidate as a cluster member. 

The remaining 8 candidates (CFHT-PL 6, 8, 14, 15, 
18, 19, 20, 22) have very low membership probabilites 
{p <20%). Of these, CFHT-PL-14 and 18 had already 
been rejected as a non-members based on the absence of 
lithium absorption in their spectrum and CFHT-PL-19, 
20, and 22 were suspected non-members based on their 
peculiar location in near-IR CMDs. The low member- 
ship probability we measure for these objects confirmed 
that they are non-members. The 3 remaining candidates 
(CFHT-PL 6, 8, and 15) deserve further discussion. 

CFHT-PL-6's proper motion is high and more than 3a 
away from the cluster's mean motion. On the one hand, 
CFHT-PL-6 does not exhibit Ha emission (Martin et al. 
2000) but this alone does not rule out membership since 
its twin, CFHT-PL-7, also lacks Ha emission but is a 
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Table 1. Proper motion of Pleiades VLM and BD candidates. The Pleiades cluster proper motion is /i Q cos(5 = 
19.15 ± 0.23 mas/yr and JlE = -45.72 ± 0.18 mas/yr. 



Name 


Utner id. 


(mas/yr) 


1 _ 

fJ.6 ± 0"Ma 
(mas/yr) 


(pL a COsS — jJ, a COsS) 
&fj, a cosS 


(fis - fig) 
°m 


CFHT-PL-1 




22.4 ± 8.5 


—48.7 ± 7.5 


0.4 


0.4 


CFHT-PL-2 




15.4 ± 9.3 


—41.1 ± 8.5 


0.4 


0.5 


flFHT-PT,^ 

Kyi. 11 _L J. U U 


HHT99 

.L J.J..I.- J Z. Z. 










CFHT-PL-4 




17.2 ± 8.6 


—48 3 ± 8 fl 


0.2 


0.3 


riFHT-PT,- 1 ! 

Kyi 11 _I_ J. U U 




18.2 ± 7.8 


—34.5 ± 7.6 


0.1 


1.5 


Kyi li J. J. u u 




5.4 it 7.4 


— 144.2 ± 7.7 


1.9 


12.8 


OFHT-PL-7 

V^J_ 11 1 J. 1 1 l 




19.2 ± 7.4 


—34.9 ± 7.8 


0.01 


1.4 


r!FHT-PT,-8 




2.21 ± 7.2 


7.31 ± 7.2 


2.35 


7.4 


CFHT-PT.-Q* 

V_ ' 1 11 1_ 1 U CI 




19.3 ± 8.4 


—31 8 + 136 


0.02 


1.0 


PFHT PT 1 fl 




98 fi 4- 8 1 


ou.u in o.o 


1 9 


1 Q 


CFHT-PT.-I 1 

Kyi 111-1 U _1_ _1_ 


Rnnup 1 


9^4 + 80 


—45.4 ± 7.8 


0.8 


0.04 


r!FHT-PT,-19 

V. 1 11 1_ 1 U l_Zi 






—43.9 ± 7.6 


2.2 


0.2 


PFHT PT.-1 ^ 


TpiHp 9 


11 Q4-7Q 

_l _i_ .o i . y 




fl Q 




CFHT-PL-14 




-0.8 ± 7.3 


-64.7 ± 7.5 


2.7 


2.5 


CFHT-PL-15 




66.5 ±8.1 


-54.1 ±8.1 


5.8 


1.0 


CFHT-PL-16 




21.7 ±7.8 


-30.8 ±7.4 


0.3 


2.0 


CFHT-PL-17 




30.8 ±8.0 


-47.8 ±7.6 


1.5 


0.3 


CFHT-PL-18 




25.0 ±7.6 


-12.9 ±7.3 


0.8 


4.5 


CFHT-PL-19 




-2.1 ±8.0 


-22.5 ± 7.9 


2.6 


2.9 


CFHT-PL-20 




-2.9 ±8.2 


16.0 ±7.9 


2.7 


7.8 


CFHT-PL-21 


Calar 3 


23.5 ±7.6 


-34.8 ±7.6 


0.6 


1.4 


CFHT-PL-22 




-20.9 ± 7.9 


-81.9 ±7.8 


5.1 


4.6 


CFHT-PL-23 




15.9 ±7.8 


-50.2 ±7.5 


0.4 


0.6 


CFHT-PL-24 


Roque 7 


24.4 ± 7.9 


-34.9 ±7.9 


0.7 


1.4 


CFHT-PL-25 




25.6 ±7.3 


-44.7 ±7.4 


0.9 


0.1 



' : CFHT-PL-3 is located close to the edge of a CCD and at the edge of the mosaic's field of view. This leads to large 
distorsions that prevent us from deriving a reliable measurement of its proper motion. 

* : The measurement error on fis is larger than average for CFHT-PL-9 due to charge transfer problems on one of the CCDs of 
the mosaic that smears the stellar profile. 



confirmed Pleiades member from its proper motion. On 
the other hand, CFHT-PL-6 lies significantly above the 
Pleiades ZAMS which led Bouvier et al. (1998) to suggest 
that it could be a nearly equal mass binary which, in turn 
might affect its short term proper motion. If it is indeed 
an equal mass binary, then the components of CFHT-PL-6 
would be substellar and therefore amenable to the lithium 
test. We also find that CFHT-PL-6 was about 0.1 mag- 
nitude brighter in the I-band in 1996 than in 1999. This 
level of photometric variablity is typical of late-M dwarfs 
(Martin et al. 1996) but might also occur in substellar ob- 
jects (Bailer- Jones & Mundt 1999). Based primarily on its 
highly discrepant proper motion and pending additional 
observations of this object, we tentatively conclude that 
CFHT-PL-6 is most likely a foreground M dwarf and not 
a cluster member. 

CFHT-PL-8 has a vanishingly small membership prob- 
ability with a proper motion close to that of background 
field stars. It had previously been considered as a proba- 
ble member based on EW(Ha) = 14. 6A, a spectral type of 
dM5.6 and IR colors consistent with membership (Martin 
et al. 2000). These properties, however, are not inconsis- 



tent with it being a low-mass field dwarf. Even the ex- 
tremely short rotational period of 0.401 days recently de- 
rived by Terndrup et al. (1999) is not unexpected for very 
low-mass dwarfs (Delfosse et al. 1998). We thus conclude 
that this candidate is not a cluster member. 

CFHT-PL-15 is puzzling in several ways. There is little 
doubt that this late-type object (Sp.T. M7) is a Pleiades 
brown dwarf since lithium has been detected in its spec- 
trum (Stauffer et al. 1998). Yet, the tangential motion we 
measure is clearly different from that of Pleiades mem- 
bers. Moreover, Stauffer et al. (1998) find that it exhibits 
unusual colors that locate it slightly below the Pleiades 
ZAMS. From the analysis of HST images, Martin et al. 
(2000) found evidence for high residuals after PSF sub- 
traction, which might indicate the presence of an unre- 
solved companion at a separation less than 0.22"and about 
3 magnitude fainter than CFHT-PL-15. It is thus conceiv- 
able that the short-term peculiar tangential motion we 
measure results from orbital motion in a binary system or 
photometric variability of one or both components. The 
latter appears more likely since, with an estimated sub- 
stellar mass, the maximum displacement of the photocen- 
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ter due to orbital motion would be less than 25 mas over 
3 years while we measure ~ 135 mas. However, CFHT- 
PL-15 does not seem to be a photometric binary from its 
location in a color-magnitude diagram. For the time be- 
ing, we thus consider CFHT-PL-15 as a highly probable 
Pleiades brown dwarf based on the presence of lithium in 
its spectrum. 

4. Discussion 

The firm identification of Pleiades VLM stars and brown 
dwarfs from their kinematics provides a clean, albeit small, 
sample of substellar objects from which a more reliable es- 
timate of the substellar mass function of the cluster can 
be derived. In addition, these are the first proper motion 
measurements obtained with an accuracy of better than 
10 mas/yr for very low mass Pleiades members and this 
allows us to start to investigate the intrinsic velocity dis- 
tribution of brown dwarfs in the cluster. These two aspects 
are briefly discussed below after the status of the objects 
contaminating the photometric sample is investigated. 

4.1. The nature of the contaminating objects 

Seven objects of the original photometric sample are found 
to be probable non members. These objects are located in 
the same region of the (R-I, I) CMD as bona fide Pleiades 
VLM stars and brown dwarfs. Their proper motions in 
Figure ^ do not seem to be uniformly distributed but 
tend to scatter over the lower left quadrant of the vector- 
point diagram. In order to investigate the status of the 
contaminating objects we used the kinematical model of 
the Galaxy developed by Robin & Creze (1986). From 
the model, we constructed a synthetic sample of field 
stars assumed to be observed in the galactic direction 
of the Pleiades cluster and covering a magnitude range 
I = 16 - 20 and a color range R - I = 1.7 - 2.7. The 
sample thus computed contains nearly 10 4 stars of which 
13 lie in the same location of the CMD as the Pleiades 
VLM and BD candidates. According to the model, these 
13 objects are late M-dwarfs (M6-M9) distributed over a 
distance between 70 and 170 pc and their proper motions 
indeed scatter preferentially over the lower left quadrant 
of the VPD with an amplitude up to 120 mas/yr. The 
same eccentric distribution of proper motions is seen for 
the whole sample of 10 4 stars and results from the ef- 
fect of galactic differential rotation in the direction of the 
Pleiades cluster. We thus conclude that the photometric 
candidates rejected as being non Pleiades member based 
on their proper motion are a mixture of foreground and 
background late-M dwarfs of the galactic disk. 

4-2. The Pleiades substellar IMF 

From accurate proper motion measurements, we have 
identified 7 out of 25 VLM and BD candidate members 



as probable field stars. In addition, the lowest mass can- 
didate of Bouvier et al.'s (1998) sample, CFHT-PL-26, 
was not included in this study but Martin et al. (2000) 
classified it as a non-member based on its spectral prop- 
erties (discrepant pseudo-continuum indices and lack of 
H a emission). The overall level of contamination of the 
photometric sample by field dwarfs is thus 8/26 = 31%, 
close to the original estimate of Bouvier et al. (1998) based 
on statistical arguments. 

Do these results modify the earlier estimate of the sub- 
stellar mass function of the Pleiades cluster? Of the 18 
photometric candidates with a mass between 0.04 and 
0.08M Q in the original sample (CFHT-PL-9 to 26), 6 
are rejected here as non-members. Thus counting 12 con- 
firmed objects in this mass range, a power law fit to 
the mass function between 0.3M Q and O.O4M yields: 
dN/dM oc M~°- 51±015 , i.e., slightly shallower than the 
earlier estimate (dN/dM oc M~ 0M0±0 - 15 ). This slight re- 
vision is not very significant considering that the errors 
are still dominated by Poisson noise from small samples. 

These results do not account for binarity. CFHT-PL- 
6 (whose membership is however in doubt, see above), 
12, and 16 appear to lie on the binary sequence of the 
cluster in a color-magnitude diagram. If these 3 objects 
are nearly equal mass cluster binaries, one finds dN/dM oc 

jy-0.65±0.15 

4--S. Kinematics of Pleiades VLM stars and BDs 

The \i a cos 6 and [1$ distributions of the Pleiades VLM and 
BD candidates are illustrated in Fig. ||. A gaussian fit to 
the distributions of the confirmed Pleiades members leads 
to dispersions of a fiaCOS s — 7.2 mas/yr and er M(5 = 8.5 
mas/yr that are similar to the expected measurement er- 
ror. Hence, we find no evidence for an intrinsic dispersion 
in the distribution of tangential velocities of the VLM and 
BD Pleiades members. This is consistent with the hypoth- 
esis that these VLM and substellar objects formed in the 
same way as stellar cluster members did, i.e., from the 
collapse of isolated, very low mass molecular cloud cores. 
Indeed, the internal dispersion of Pleiades stars is ~ 1 
mas/yr (Van Leeuwen, 1980). 

However, owing to the difficulty of triggering the col- 
lapse of molecular cores whose mass is much lower than the 
typical Jeans mass (~ O.7M , Clarke et al. 2000), alterna- 
tive models have been proposed for the formation of brown 
dwarfs. Burkert et al. (1997) have shown that collapsing 
molecular cores are prone to multiple fragmentation that 
eventually leads to the formation of small-N protoclusters, 
including a number of very low mass fragments. Following 
dynamical interactions within the protocluster, the least 
massive fragments are preferentially ejected with typical 
velocities of order of a few km/s (Burkert, priv. comm.) 
and may thus become isolated brown dwarfs (Klessen & 
Burkert 2000). Alternatively, Lin et al. (1998) proposed 
that isolated brown dwarfs may form as a result of an 
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Table 2. Candidates membership criteria. The EW(Li) are from Stauffer et al. (1998) and the others tests have been 
performed by Martin et al. (2000). In column 8, we give the membership probability of the candidates based on their 
proper motion. Our final assessment regarding membership is indicated column 9. 



Name 


other id. 


EW(Li) 
(A) 


H a 


Vrad 


I-K 


Martin et al. 


Membership 
prob. (%) 


Conclusion 


CFHT-PL-1 






yes 




yes 


yes 


98.6 


yes 


CFHT-PL-2 






yes 




yes 


yes 


98.9 


yes 


CFHT-PL-3 


HHJ22 








yes 


yes 




yes 
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Fig. 2. Proper motion distributions for program stars. 
The hatched histogram refers to candidates confirmed as 
being Pleiades members and the solid curve is a gaussian 
fit to their proper motion distribution. 



encounter between protostars with massive disks. The en- 
counter leads to the formation of an unbound tidal tail 
which contains part of the initial disk's mass and may 
later condensed to form an isolated substellar object. Here 
again, the typical ejection velocities are of order of a few 
km/s. 

At a distance of 118 pc, a tangential velocity of 1 km/s 
corresponds to a tangential motion of about 2 mas/yr. 
Proto brown bwarfs ejected with this velocity in the early 



stages of cluster formation, about 120 Myr ago, would 
have now drifted by several tens of degrees away from 
the cluster. If these BD formation models are correct, the 
sample we have observed close to the cluster's center must 
represent a tiny fraction of the primordial BDs, namely 
those which populated the low tail of the distribution of 
ejection velocities. In such a case, the slope of the substel- 
lar mass function of the Pleiades cluster might currently 
be largely underestimated. 

5. Conclusion 

A sample of very low-mass stars and brown dwarfs have 
been firmly identified in the Pleiades cluster on the basis 
of their proper motion. With 12 confirmed objects with 
a mass less than 0.08-M© distributed over an area of 2.5 
square degrees close to the cluster's center, the Pleiades 
mass function is estimated to be dN/dM oc M°- 51±al5 in 
the mass domain ranging from 0.04 to 0.3M Q . Taking in- 
dividual components of suspected substellar binaries into 
account leads to dN/dM oc M (35±015 . The main source 
of uncertainty on these estimates now lies in the unknown 
radial distribution of brown dwarfs relative to stars in 
the cluster. As a group, the identified brown dwarfs ex- 
hibit no intrinsic velocity dispersion. This suggests that 
they have formed from the collapse of isolated low mass 
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molecular cores. However, we cannot rule out that this 
sample of brown dwarfs represent only a tiny fraction of 
the brown dwarf population originally formed in the clus- 
ter and which might have escaped since then either due 
to dynamical ejection at the time of their formation or 
during the subsequent dynamical evolution of the cluster. 
Additional studies covering a much larger fraction of the 
cluster's area are needed to settle this issue. 
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